Thermal spin flips in atom chips 
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We derive a general expression for the spin-flip rate of an atom trapped near an arbitrary dielectric 
body and we apply this theory to the case of a 2-layer cylindrical metal wire. The spin flip lifetimes 
we calculate are compared with those expected for an atom near a metallic slab and with those 
measured by Jones et al. above a 2-layer wire [M.P.A. Jones, C.J. Vale, D. Sahagun, B.V. Hall, and 
E.A. Hinds, Phys. Rev. Lett. 91, 080401 (2003)]. We investigate how the lifetime depends on the 
skin depth of the material and on the scaling of the dimensions. This leads us to some conclusions 
about the design of integrated circuits for manipulating ultra-cold atoms (atom chips). 
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I. INTRODUCTION 

Microscopic traps provide a powerful tool for the con- 
trol and manipulation of Bose-Einstein condensates over 
micrometer distances. Microstructured surfaces, known 
as atom chips, are particularly interesting for this pur- 
pose since they can be tailored to provide a variety of 
trapping geometries T] and promise well-controlled quan- 
tum state manipulations of neutral atoms in integrated 
and scalable microtrap arrays. Ultimately there is the 
possibility of controlling the quantum coherences within 
arrays of individual atoms for use in quantum informa- 
tion processing Q ■ This technology is attractive because 
it appears robust and scalable and because trapped neu- 
tral atoms can have long coherence times. 

However, atoms in these traps are held close to the 
micro-structured material surfaces, which are typically 
at room temperature. Thermal fluctuations give rise to 
Johnson noise currents in the material Q . Such currents 
are normally observed as a noise voltage across a resis- 
tor, but they also cause the electromagnetic field near 
a conducting solid to fluctuate with a broad noise spec- 
trum. For atoms trapped close to the surface of a con- 
ductor these fluctuating fields can be strong enough to 
drive rf magnetic dipole transitions that flip the atomic 
spin. If the atom is in a magnetic trap where only low- 
field-seeking Zeeman sublevels are confined, the spin flips 
lead to atom loss. This is known experimentally [i|,|5j as 
well as theoretically @. The loss rate increases strongly 
as the atoms approach the metallic surface of an atom 
chip. For a given desired lifetime, this restricts how close 
the trapped atoms can be brought to the surface, which 
in turn determines the period of the smallest trapping 
structures that can be imposed on the atom by the chip. 

The paper is organized as follows: In Section^ we in- 
troduce the basic equations and discuss the quantization 
of an electromagnetic field in the presence of a dispersing 
and absorbing dielectric body. Then, in Section IlIII we 



derive a general expression for the spontaneous and ther- 
mal spin-flip rates of an atom due to the coupling of its 
magnetic moment to the magnetic field. This derivation 
is based on the Zeeman Hamiltonian of the system and 
the corresponding Hciscnbcrg equations of motion. We 
show that the spin-flip rate is determined by the dyadic 
Green tensor of the classical, phenomenological Maxwell 
equations. In Section IdTI we present the scattering Green 
tensor for a 2-layer cylindrical body surrounded by an un- 
bounded homogeneous medium, with details being given 
in Appendix ^ Then, in Section ^ we use this Green 
tensor to obtain an explicit analytical expression for the 
total spin-flip rate of an atom above a 2-layer wire. Some 
numerical results are presented and discussed in Section 
IVII The numerical results are compared with the cor- 
responding results for a slab and with the experimental 
measurements presented by Jones et al. in Ref. [4]. Our 
conclusions are given in Section IVIII 



II. BASIC EQUATIONS AND QUANTIZATION 

In classical electrodynamics, dielectric matter is com- 
monly described in terms of a phenomenologically intro- 
duced dielectric susceptibility. Let us consider a classical 
electromagnetic field, described by the phenomenologi- 
cal Maxwell's equations, without external sources. We 
restrict our attention to isotropic but arbitrarily inhomo- 
geneous non-magnetic media, and assume that the polar- 
ization responds linearly and locally to the electric field. 
A linear response formalism similar to that presented be- 
low can also be found in Refs. 

The most general relation between the matter polar- 
ization and the electric field consistent with causality and 
the fluctuation-dissipation theorem is [j| 



P(r,t)=e J dr x{ r , t) E(r, t — r)+ P»(r, t) , (1) 



* Electronic address: |s . scheeKBimperial . ac . uk| 



where x(r, i) is the linear susceptibility. The inclusion 
of the noise polarization P^(r,t) is necessary to fulfil 
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the fluctuation-dissipation theorem. It is this fluctuating 
part of the polarization that is unavoidably connected 
with the loss in the medium. Converting the displace- 
ment field D(r, t) = £oE(r, t) + P(r, t) into Fourier space 
using Eq. (JIJ, we obtain 

D(r, lu) = e £(r, w)E(r, u) + Pjv(r, w) , (2) 

where e(r,o;) is the complex permittivity and s(r, uj) — 1 
is the temporal Fourier transform of %(r, f). The real 
part of the permittivity (sr, responsible for dispersion) 
and the imaginary part (ej, responsible for absorption) 
are related to each other by the Kramers-Kronig relation. 

Using Maxwell's equations in Fourier space, we find 
that E(r, uj) satisfies the Hclmholtz equation 

uj 2 

V x V x E(r,w) 5-e(r, u;)E(r, uj) = uj 2 n Q P N (r,uj) , 

c z 

(3) 

with the solution 



E(r, uj) = wVo J d 3 r' G(r, r', w) • P^(r', w) , 



(4) 



where the Green tensor G(r, r',u>) is a second-rank tensor 
determined by the partial differential equation 

UJ 2 

V x V x G(r, r',u)- — e(r, u)G(r, r',oj) = S(r -r')U, 

where U is the unit dyad. Together with the boundary 
condition at infinity, this equation has a unique solution. 
In accordance with Maxwell's equations the correspond- 
ing solution for the magnetic field in Fourier space is 
B(r,w) = (iuj)- 1 V x E(r,w) . 

As we have seen, the noise polarization P n (r, t) plays 
a fundamental role in determining the electric field. The 
form of Pjv (r, t) follows from the fluctuation-dissipation 
theorem, which states that the fluctuations of the macro- 
scopic polarization are given by the imaginary part of the 
response function [here £j(r, uj)]. If we pull out a factor 
and define the dynamical variables f (r, uj) as the funda- 
mental (5-correlated random process, we find that we can 
write the noise polarization as Q 



P \ ( i ^ ! = / \/ ^ £/(r, uj) f (r, uj) . 



(6) 



Upon quantization, we replace the classical fields f (r, uj) 
by the operator- valued bosonic fields f(r, w) which we 
associate with the elementary excitations of the system 
composed of the electromagnetic field and the absorbing 
dielectric matter. They satisfy the equal-time commuta- 
tion relations [/j(r, uj), /J(r', uj')] = Sij5(r — r')8(u> — uj'). 

The magnetic-field operator in the Schrodinger picture 
can now be obtained as 

B(r) = B(+)(r) +B(-)(r), BH( r ) = [B(+)(r| (7) 

where 



oo 

B (+) (r) = J du>B(r,uj) , 



is its positive- frequency part. In this way, the electro- 
magnetic field is expressed in terms of the classical Green 
tensor satisfying the Helmholtz equation (JHJ and the con- 
tinuum of the fundamental bosonic field variables f (r, uj). 
All the information about the dielectric matter is con- 
tained, via the permittivity e(r,uj), in the Green tensor 
of the classical problem. 

We close this Section by mentioning two important 
properties of the Green tensor. It can be shown that the 
(Onsager) reciprocity relation G(r, r' ,u>) = G T (r',r,uj) 
holds [l0||. Additionally, another useful property is the 
integral relation 

Jd 3 r' ^ £l (v' .u)G kl (v,v' ,uj)G* nl {T A y ,uj) (9) 

= ImGt„(r,r A ,w) , 

which we will use later in this paper. Both relations es- 
sentially follow from linear response theory, with Eq. @ 
being equivalent to the fluctuation-dissipation theorem 
[Tl|. It should be noted that we assume the dielectric 
permittivity to possess at least an infinitesimal imagi- 
nary part everywhere to avoid surface contributions in 
Eq. ©. 



III. DERIVATION OF THE SPONTANEOUS 
AND THERMAL SPIN FLIP RATES 

The Hamiltonian of the combined system of electro- 
magnetic field and absorbing matter, from which the 
(quantized) phenomenological Maxwell's equations can 
be derived, can be written in terms of the basic field op- 
erators f(r, uj) in the diagonal form 

oo 

H = J d 3 r J dufou) ft(r,u) ■ f(r,w) + ^ ftwa| Q , 
o a=i -f 

(10) 

which leads in the Heisenberg picture to the (quasi-free) 
time evolution f(r, uj) — > f (r, ui)e~ %ult . Here we have also 
included an atom through the operators £ Q = |a)(a| and 
the energy huj a of the atomic state \a) (a = i, /). 

The interaction of the atom at position ta with a mag- 
netic field B(r) is described by the Zeeman Hamilto- 
nian Hz = —ft ■ B(i-a), where fi = fi\i)(f\ + h.c. is the 
magnetic moment operator associated with the transition 
\i) — ► I/). The magnetic moment vector is 



i\ [i B gsS + .9lL - gi 1 \f) , 



(11) 



(8) 



where /i_b is the Bohr magneton, S is the electronic spin 
operator, L is the orbital angular momentum operator, 
I is the nuclear spin operator and gs ~ 2, g^ and gi are 
the corresponding ^-factors. We restrict our attention 
to L = 0, which corresponds to the ground state of an 
alkali atom, and we neglect the small nuclear magnetic 
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moment in comparison with the Bohr magneton. In the 
rotating-wave approximation, we can then write the Zee- 
man Hamiltonian as 



-VB9S 



(f\S q \i)^B^(r A )+h. C .] , (12) 



where the atomic raising (lowering) operator = 
\i)(f\ = (C*- - - 1 )^] satisfies the commutation relation 

= ±i (± \ with l z = - |/) (/|). Repeated 

indices q indicate a sum over spatial vector components. 

Using the Hamiltonian JHJ, the Heisenberg equation 
of motion for the atomic quantity £ z (t) is given by 



VB9S 

ih 



(f\S q \i)i^B(+\r A )+h.c.. (13) 



Furthermore, the Heisenberg equation of motion for the 
bosonic field operator is 



/»(r,w,f) 



(14) 



{i\Sq\f)i { )e q P j9p-oy/si(r, a) GUr A , r, uf) 



where e qp j is the Levi-Civita symbol and dj = d/dxj. 
This equation can now be formally integrated to yield 



/ < (r,w J t)=/ i , fre e(r,w J t)+ dTe-^^^-Hr) (15) 



WBgs 



i\S q \f) £qp 3 d p -^y / e I (r,u;)G* l (r A , r,u), 



where /i,f re e( r > t) denotes the freely evolving basic-field 
operators. The lowering operator ^'(t) in Eq. (fT5|l can 
be found by solving its Heisenberg equation of motion. In 
the Markov approximation, this solution can be reduced 
to its slowly varying part £(~^(t) e 4Wi/ (* _T ' in Eq. I|15|l so 
that the time integral can be approximated by 

/i(r jW ,t) = jUe^t) + ^gE= 



Xe W j d p -^y/ei(r, uj) G* t (r A , r, w) C(w</ - w) , (16) 

where £(x) = 7t<5(:e) + iVx^ 1 (V denotes the principal 
value) and u>if = u>i — u>f is the transition frequency cor- 
responding to the flip |i) — > |/) in the atom's internal 
state. Substituting this formal solution into the expres- 
sion for the magnetic field, we obtain 

B(+)(r A ,u>,t) = B%l ee (r A ,u>,t) 

+ B S ° (i|5 , p|/)e 9 jfee P mn9 : ,i9 m ^")(i)C(w i/ -w) 

7T 

x d 3 r — e I {r,u;)Gki(r A ,r,uj)G* ll {r A ,r,uj). (17) 

The spatial integral can be evaluated using the integral 
relation Eq. 10 yielding lmGk n (^A, ?a, w). Therefore, 



Eq. |JT7|| becomes 



B(+>(r^w,t) = flW (r A ,w,t) 



Im[ V x V x G(r^,rA,w) 



(18) 



Performing the ^-integration and inserting into Eq. l|13fl . 
we obtain 



iSlu ' 



+ 



(f\Sq\i)£ W B^l e (v A ) + h.c. 



(19) 



where the spontaneous spin- flip rate T B = T B (r A ,u>if) 
arises from the 5 function (the real part of the £ function) 
and is given by 



f b = ^He^ {fl s qm s p \ f) 

x Im[ V x V x G(r A ,r A ,u)i f ) } 



qp 7 



(20) 



and where the term Su arises from the principal-value 
integral (the imaginary part of the £ function) and is 
identified as the radiative frequency shift. Furthermore, 
the shifted frequency is given by cDy = toij + 5uj. In what 
follows, the transition frequency is always taken to be the 
shifted frequency Ui / that one measures in an experiment 
and not the bare frequency u>if = uj. For simplicity, we 
omit the tilde in all subsequent formulas. Note that the 
same result for T B is obtained when using an appropri- 
ately derived master equation as done in Ref. |6j. 

We assume that the dielectric body is in thermal equi- 
librium with its surroundings. The magnetic field is then 
in a thermal state with a temperature T, equal to the 
temperature of the dielectric body. The total flip rate 
for the atom is therefore given by r^ tal = r B (nth + 1), 
where the mean thermal occupation number is 



1 



nth 



1 



(21) 



and /cb is Boltzmann's constant. At zero temperature, 
i.e. n t h = 0, the relaxation dynamics is entirely due to 
the spontaneous flip rate T B . For large T on the other 
hand, n t h ~ i^- ~> 1 and the spin flip rate is predomi- 
nantly induced by thermal fluctuations. 

In the experiment of Ref. Q 87 Rb atoms are initially 
pumped into the trapped state \F,m) = |2,2). Thermal 
fluctuations of the magnetic field then cause the atoms 
to evolve into hyperfine sublevels with lower tuf- Upon 
making a transition to the mp = 1 state, the atoms are 
more weakly trapped and are largely lost from the re- 
gion of observation, causing the measured atom number 
to decay with a rate . Here we are introducing the no- 
tation T B l . nlf for the total spin-flip rate associated with 
the transition |2, rrij) — > |2,m/). 
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IV. THE DYADIC GREEN TENSOR 

The geometry we are considering in this paper is a 2- 
layer cylinder surrounded by an unbounded homogeneous 
medium (see Fig. This corresponds to the experimen- 
tal geometry in Ref. Q . Because the Helmholtz equation 



layer 3 
(vacuum) 




layer 2 



layer 1 



atom 




FIG. 1: The geometry we are considering is a 2-layer cylin- 
der surrounded by an unbounded homogeneous medium. The 
outer region is labelled layer 3 (vacuum), the coating is layer 2 
and the cylinder core is layer 1. The distance from the surface 
of the outermost layer to the atom is r. 

is linear, the associated Green tensor can be written as a 
sum, 

G(r, r', u) = G°(r, r', u) + G wirc (r, r', u) , (22) 

where G°(r,r',w) represents the contribution from the 
vacuum and G wlrc (r, r', uj) describes the part due to the 
wire. When the atom is located in layer 3, the scattering 
contribution is 

EI 



G wirc (r,r» = — 



where 



dh 



ra=0 



P^K n (h), (23) 



W£ 3 



(24) 



N^(h)M^(-h)-NW(h)M 



-h) 



+M^(h)N'W(-h) M<$(h)NZ\-h) 



Rf{h) 



M^(h)M£\-h) + M$>(h)M'M(-h) 



For simplicity, we have omitted the tensor product sym- 
bol ® between the even and odd cylindrical vector func- 
tions defined by Me„(/i) = V x [ipe n (h)z] and Ne n (h) = 
V x V x [ipe n (h)z]/k 3 . The scalar eigenfunctions ]ju n {h) 
satisfy the homogeneous scalar wave equation [12j. It 
follows from these definitions that 



Explicitly, 



V x M %n (h) = k 3 N. n (h) , 

V x N %n (h) = k 3 M. n (h). 



dZ n (i] 3 p) cos , . 

in : . {n<p)e p 

dp sin v ^' p 



Tih-Z n (rj s p) ^ (n<?!))e + T)lZ n (r] 3 p) (ncf>)e 2 



(25) 
(26) 



(27) 



Ahz 



MeJh) 



n . . sm , 
T -Z n (r) 3 p) _ c (n^)e p 



cos 



dZ n (rj 3 p) cos , 

. (nffl)e^ 

dp sm v ^ 9 



Ahz 



(28) 



The primes in Eq. I|24|) indicate the spherical coordinates 
(p' ,(/)', z'). The superscript (1) indicates that Z n should 

be replaced by the Hankel function of first kind H„ . 
Otherwise, Z n is the Bessel function of first kind J n . The 
propagation constant in the p direction is rj^ = k p — /i 2 , 
where k p is the wave number of the pth layer. The per- 
mittivity of the pth layer is denoted by e p . The scatter- 
ing reflection coefficients R„(h) are given in Appendix lAl 
(k, 1 = 1,2). 

The double curl of the Green tensor in Eq. I|23(l can be 
written 



V x V' x G wir ' 



n=0 



(-^n)xx {J-Th)xy {Jn)xz 
{ln)yx [Injyy {^n)yz 
{Ifijzx \^n)zy {Jn) zz 



(29) 



where 

,00 1 

I„ = I„(r,r',u/) = / dh -j V X V X n n {h) . (30) 

Note that the curls are computed by replacing N« n (ft) 
by Me„(/i) and vice versa, according to Eqs. 
and Also note that the integration variable h is 

the wave number in the z-direction (see Fig. ^) . From 
the symmetry of the integrand, it is easy to show that 
' ) zx = (Pr)y Z = (Ilr)zy = 0, where 
(I n (r,r',cj))ij = x,y,z). Note 
that the (Onsager) reciprocity relation as mentioned in 



{I n )xz — {I n • "' 



Section CU implies that (4 in %- = (4 im )j 8 . 



V. THE SPIN FLIP RATE OUTSIDE A 2-LAYER 
WIRE 

The spin-flip rate in free space is readily derived 
from Eq. 1(20)1 since Im[V x Vx G°(ta, ta, ^if)]qp = 



5 



(fc|/67r) 8 qp , where k 3 = uj/c is the free space wave num- 
ber corresponding to the atomic transition. We use the 
notation k 3 here because in our discussion of the cylin- 
drical wire, the third layer is a vacuum. Hence 



( mb.9s) 2 

37T H 



1,3 q 2 
K 3°if ' 



(31) 



where we have introduced the angular factor = 

mS x \f)\ 2 +|<i|4|/)| 2 = 2Ki|4|/)| 2 - We do not have a 
term containing S z since we are interested here in a spin 
flip, which by definition changes mp. We have moreover 
used the fact that the two transverse matrix elements are 
equal in absolute value as a result of symmetry. For the 
87 Rb ground state transition \F tuf) = |2,2) — > |2, 1), 
the angular factor S^ 2 = 1/8 [lS| . 

In order to find the contribution of the wire to the 
spontaneous spin-flip rate, we use Eqs. I|2UI) and (|29H . 
The quantization axis is taken to be along the z direction, 
corresponding to the direction of the bias field that the 
trapped atoms experience in the experiment. We obtain 



We see from Eq. ()22l) that the total spin-flip rate is 
equal to the sum of the free space contribution and the 
scattering contribution. The total spin-flip rate for the 
rate-limiting transition |2,2) — > |2, 1) is therefore 



-pwirc 
1 21 , 



fab + 1) 



(35) 



VI. NUMERICAL RESULTS 

In the experiment of Ref. |3| , cold atoms are held in a 
microscopic trap near a current-carrying wire assumed 
to be at room temperature. The lifetime for atoms to 
remain in the microtrap is measured over a range of 
distances down to 27 /im from the surface of the wire. 
The wire consists of a central copper core with radius 
rii = 185 ^m and a 55 /im thick aluminium layer, i.e. 
a 2 = 240 /im. Using Eq. 13411 . the resistivities p\ = 
1.6 ■ 1CT 8 Om for Cu and p 2 = 2.7 • 1CT 8 17m for Al give 
skin depths of Si = 85 /im for Cu and 5 2 = 110 /im for Al 
at frequency / = ui/2tt = 560 kHz. 



THwire T-i \ ^ / 



(2-<W)Re (?„ lm ) ra + (/r )vv -(32) 



fin) 



Here we have once again used the facts that K^S^I/)! 2 = 
I (* I A/ 1/) 1 2 and that (i 1,5*1/) = 0. The dimensionless in- 
tegrals (7£ m )y = {I^ij/kl are given by 



or) 



dq x 

1 m 



R^(q)+ q 'Rtf( q )][(H n3 ) 



k^(a 2 +r) 2 



(33) 



+ (mK 3 ) 2 } 



k 3 (a 2 + r) 



since p 



a 2 



We have used the simplified notation 
Z n (jjpk 3 p) and the primes in Eq. Ij33(l denote the 
derivative with respect to the full argument of the rele- 
vant function, e.g. Z' np = dZ n {rj p k 3 p) j 1 d{r\ p k 3 p). We have 
also chosen to write the permittivity of the pth layer rel- 
ative to the outermost layer, i.e. e 
number for layer p is then given by k p = k 3 e p 
the dimensionless propagation constant rj p = i] p /k 3 in 

2 where 



p — £3£p Cl . The wave 
2 _ f.2^01 and 



the p direction can be written as rip = W e r p cl — q 

q = h/k 3 is the dimensionless integration variable. 

The skin depth is the characteristic length scale on 
which an electromagnetic wave is damped within a con- 
ducting medium. It is given by S p = ^J2e^p p uj /fco (see 
e.g. Ref. where p p is the resistivity of layer p. Since 

the spin-flip frequency is very much lower than the reso- 
nance frequencies of the material in the wire, the relative 
permittivity is related to the skin depth by [l3| 



_rel 



EoPpUJ 



n o p 



(34) 




40 60 

r (/im) 



100 



FIG. 2: Lifetime r of the trapped atom as a function of 
the atom-surface distance r. Dotted curve: calculated spin- 
flip lifetime near a 2-layer wire at 300 K with the parameters 
/ = 560 kHz, a± = 185 pm, a-2 = 240 pm, Si — 85 pm, and 
^2 = 110 pm. Solid curve: The same but at 380 K. Dot-dashed 
curve: calculated lifetime near a thick Al slab at 380 K with 
S = 110 pm (using Eq. (35) of Ref. @). Crosses: measured 
lifetimes of Ref. y| . 

The dotted line in Fig. shows the lifetime r = 1/rfj 
that we have calculated assuming a temperature of 300 K, 
together with the measured lifetimes (crosses). We see 
that the experimental results are close to the theory, in- 
dicating that the thermal spin flip mechanism is the pri- 
mary cause of atom loss in the experiment. Neverthe- 
less there is also a clear systematic discrepancy, with the 
measured lifetimes being 20 — 30% shorter than expected. 
We find excellent agreement when the temperature in our 
theory is increased to 380 K, as shown by the solid curve 
in Fig. |2 We have re-examined the conditions under 
which the experiment was run and consider it most likely 



6 



that the wire temperature was indeed ~ 380 K, rather 
than the 300 K previously assumed. Such a temperature 
rise would be consistent with known power dissipation 
and with reasonable assumptions about the heat flow. In 
effect, the thermally driven spin flips have allowed us to 
measure the temperature of the wire! 

The theory for the decay rate of an atom above a plane, 
thick slab is already known Applying this theory to 
an Al slab with skin depth 8 = 110 /mi and temperature 
380 K, we obtain the result shown dot-dashed in Fig. [21 
This curve lies below that for the wire, simply reflecting 
the fact that the slab contains a larger volume of fluc- 
tuating polarization than the wire. Naturally, the two 
380 K curves converge at sufficiently small atom-surface 
distances (r <C 62, Si, 0.2), and in that range they vary 
linearly with distance [(| . 

100 r . . 1 i 




1 1 1 1 1 1 

1 10 100 1000 

82 (/im) 

FIG. 3: Lifetime r of the trapped atom as a function of the 
skin depth $2 of the outer layer. The atom-surface distance 
is fixed at r = 50 /im. The other parameters are: / = 560 
kHz, T = 300 K, ai = 185 /mi, a 2 = 240 /mi, and 5i = 85 /mi. 
The straight dashed line represents the large 82 limit. The 
numerical value for this limit is 52 s. 

The lifetime for the atom to remain in the trap exhibits 
a minimum with respect to variation of the skin depth, 
as illustrated in Fig-EI where the skin depth of the wire 
core 81 is fixed at 85/im but the skin depth of the outer 
layer 82 is varied. Below the minimum at 82 — 20 fim, 
a decrease of skin depth leads to an increase of lifetime 
in proportion to 5^ ■ This happens despite a growth in 
the polarization noise [see Eqs. 10 and (|34|l ] because 
the region generating the noise is becoming thinner. In 
the small 82 limit the outer layer approaches a perfect 
conductor, the core wire does not play any role, and the 
lifetime becomes exceedingly long. By contrast, when the 
skin depth increases above 20 /mi, the reduction in polar- 
ization noise is more influential than the growth of the 
source volume. In this region it is the worse conductor 
that gives the longer lifetime. At large 82, the outer layer 
of the wire approaches the free space limit, and the life- 
time is entirely determined by the skin depth and radius 
of the core. From a practical viewpoint it would normally 



be desirable to avoid the minimum of the lifetime curve. 
This means avoiding surface materials whose skin depth 
at the spin flip frequency is comparable with the atom- 
surface distance. This is a generic result. For example, 
at height z above a slab, one obtains the shortest lifetime 
when the skin depth is z/3 1 ^ 3 , as is readily derived from 
equation (23) of Ref. Q. 




0.1 1 1 — 1 — 1 — 1 

1 10 100 1000 

r (/im) 

FIG. 4: Lifetime r as a function of atom-surface distance r, 
with r, ai and 0,2 scaling together according to a 2 = 5r and 
ai = (185/240)a2- The other parameters are: / = 560 kHz, 
T = 300 K, 81 = 85 /mi, and 82 = 110 /im. 

From the same perspective of cold atoms trapped 
above small integrated circuits (atom chips) it is also 
interesting to see how the lifetime is altered when the 
dimensions 01, a 2 of the wire are varied or the atom- 
surface distance r changes. For example, let us scale 
all three lengths together, such that r = 0,2/5 and 
ai = (185/240)a2, while the skin depths are fixed. The 
result of such a scaling is illustrated in Fig. 0] When the 
atom-surface distance is large compared with the skin 
depth, i.e. r = 0,2/5 S> 62 ~ 100 /im, the spin-flip life- 
time scales as ~ r 4 . This has the same exponent as the z 4 
scaling of lifetime that applies at distance z from a slab 
in the range where z 3> 8 |fj . The correspondence seems 
natural to us since the wire is essentially a curved slab 
when the skin depth is small. For a given ratio of atom- 
surface distance to wire size, the two lifetimes should 
therefore be related by a constant geometrical factor, re- 
sulting in the same distance scaling. At the opposite 
extreme, where z <C 6, the slab result is r oc z. By con- 
trast, we see in Fig. 0] that the lifetime outside the wire 
approaches a constant when 02 = 5r -C 82- This differ- 
ence occurs because the thickness of the source region is 
not the skin depth, but rather the diameter of the wire, 
which we are scaling linearly with r. In a similar way, 
it is possible to lengthen the lifetime of an atom above 
a slab by reducing the thickness of the slab to less than 
the skin depth 0]. 

As a second example of scaling, we change the diame- 
ter of the wire, keeping a\ = (185/240)et2 but fixing the 
distance from the surface at 50 /tm. Once again the skin 
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depths are fixed. The resulting variation in the lifetime 
of the atom with wire size is shown in Fig. [5J At large 
wire diameter, the lifetime approaches 8.2 s, which is of 
course the same as the lifetime 50 /im above a slab with 
110/im skin depth. By contrast, when the wire size is 
small, i.e. a 2 <C r, 62, the decreasing volume of mate- 
rial leads to a a^ 3 scaling of the lifetime. In the limit 
a,2 — > 0, r™ rc vanishes and the lifetime for the atoms to 
remain in the trap is just the free-space blackbody rate 
given by the first term in Eq. (|35|) . For / = 560 kHz 
and T = 300 K, this free space lifetime is an astonishing 
- 10 18 s (see also 0). This figure emphasizes the very 
low strength of the electromagnetic field fluctuations in 
free space compared with those near a dielectric medium 
due to the surface modes. 

10000 f 1 , , 1 , , j 



the generic result that there is a minimum in the life- 
time when the skin depth is comparable with the atom- 
surface distance. When the dimensions of the wire and 
the atom-surface distance r are varied together, the life- 
time scales as r 4 at large r, following the same scaling law 
as a corresponding plane, thick slab, whereas the lifetime 
approaches a constant at small r. If instead we fix the 
atom-surface distance and vary only the dimensions of 
the wire, the lifetime scales as r -3 when the wire is small, 
leaving only the very weak free-space decay rate in the 
limit of a vanishing wire diameter. The main conclusion 
for atom chip design is that one should avoid a material 
whose skin depth at the spin flip transition frequency is 
comparable with the atom-surface distance. The lifetime 
can also be improved by making sure that metal films on 
the surface are thinner than the skin depth. 
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FIG. 5: Lifetime r as a function of outer wire radius 02 
with the atom-surface distance fixed at r = 50/im. The in- 
ner radius is scaled according to a\ = (185/240)a2. Other 
parameters are: / = 560kHz, T = 300K, 81 — 85/im, and 
82 — llOpim. Dotted line: the large ai limit. 



VII. CONCLUSIONS 

In this paper we have derived the magnetic spin flip 
rate for an atom close to an absorbing dielectric body. 
The rate is given in terms of a dyadic Green tensor, allow- 
ing the expression to be applied in principle to a dielectric 
body of any shape. We derive an explicit expression for 
the spin-flip rate of an atom outside a 2-layer cylindrical 
wire, as used in the experiment of Jones et al. p|. We 
compare our numerical results with their measurements 
and we find lifetimes marginally longer than those mea- 
sured in the experiment. The most likely explanation for 
this discrepancy is that the wire was hotter than pre- 
viously thought. We also compare the cylindrical case 
with that of a slab and show that the spin-flip lifetime 
is systematically longer above a cylinder, as one would 
expect. 

We have investigated how the lifetime of the atoms 
depends on the skin depth of the material. We find 



APPENDIX A: THE SCATTERING 
REFLECTION COEFFICIENTS 

The scattering reflection coefficients for a cylindrical 
geometry can be computed for any number of layers (see 
e.g. Refs. [H[ll E E3). In this Appendix we present 
the explicit expressions for the scattering reflection coef- 
ficients corresponding to our 3-layer cylindrical geometry. 
To find these reflection coefficients we have used the it- 
eration tensor equations in Ref. |12|. These iteration 
equations are given for arbitrary complex permittivity 
e p and arbitrary complex permeability fx p . Therefore, 
the reflection coefficients presented in this Appendix ap- 
ply to arbitrary e p and arbitrary n P . However, we stress 
that the theory presented in the main body of this pa- 
per is particular to non-magnetic media; we assume that 
Mp = Mo m all the layers p. 

The reflection coefficients are given as follows: 
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where 
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Moreover, we have 
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The function a 2 i, a 22 , and 6 2 i, & 22 are obtained from ai 2 , 
On, an d ^I2j &ii, respectively, by replacing fi p <-> — e p . 



In the last step in Eqs. (|A11|) and (|A13|I we have used 
the Wronskian determinant between Bessel and Hankcl 
functions. Finally, we have 



tn21 — 
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Whenever the combination Z 2 and Zi is involved in the 
superscript, the radius ai is implicit in the cylindrical 
functions. For example, in Eqs. (|A18(1 - IIA22(1 we have 



Z n i = Z n (rna\) , Z n2 = Z n {j] 2 ai) , 
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H 3 J 2 „ J 3 J 2 n H 3 J 2 n H 3 H 2 _^3^2 UH 3 J : 



^nu 3 , Linu 3 , une 3 , ^ne 3 , 



u ne 3 , , 

and b^ aH2 arc defined analogously, where we understand 
that the radius a 2 is implicit in all those functions. Of 
course, for the special case fi p = 1 for all layers p, these 
reflection coefficients simplify. 

The reflection coefficients R^(h) and R^(h) can be 
obtained from R n 2 (h) and R^(h), respectively, by re- 
placing /i p <-> — £ p . Note that the scattering coefficients 
R}i(h) as well as R 22 (h) are dimensionless. However, the 
coefficients R} 2 {h) and R^(h) are not, but the particular 
combinations R n 2 {h)(—ujEz/kz) — R 2 } (h)(uj fi 3 / 'fc 3 ) are. 



[1] E.A. Hinds and LA. Hughes, J. Phys. D: Appl. Phys. 32, 

R119 (1999); R. Folman, P. Kriiger, J. Schmiedmayer, 

J. Denschlag, and C. Henkel, Adv. At. Mol. Opt. Phys. 

48, 263 (2002). 
[2] T. Calarco, E.A. Hinds, D. Jaksch, J. Schmiedmayer, 

J.I. Cirac, and P. Zoller, Phys. Rev. A. 61, 022304 (2000). 
[3] J.B. Johnson, Phys. Rev. 32, 97 (1928); H. Nyquist, 

Phys. Rev. 32, 110 (1928). 
[4] M.P.A. Jones, C.J. Vale, D. Sahagun, B.V. Hall, and 

E.A. Hinds, Phys. Rev. Lett. 91, 080401 (2003). 
[5] D.M. Harber, J.M. McGuirk, J.M. Obrecht, and 

E.A. Cornell, J. Low. Temp. Phys. 133, 229-238 (2003). 
[6] C. Henkel, S. Potting and M. Wilkens, Appl. Phys. B 69, 

379 (1999); C. Henkel and M. Wilkens, Europhys. Lett. 



47, 414 (1999). 

[7] L. Knoll, S. Scheel, and D.-G. Welsch, QED in dispers- 
ing and absorbing media, in Coherence and Statistics of 
Photons and Atoms, ed. J. Pefina (Wiley, New York, 
2001); S. Scheel, L. Knoll, D.-G. Welsch, and S.M. Bar- 
nett, Phys. Rev. A 60, 1590 (1999); S. Scheel, L. Knoll 
and D.-G. Welsch, Phys. Rev. A 60, 4094 (1999); Ho 
Trung Dung, L. Knoll and D.-G. Welsch, Phys. Rev. A 
62, 053804 (2000). 

[8] G.S. Agarwal, Phys. Rev. A 11, 230 (1975); J.M. Wylie 
and J.E. Sipe, Phys. Rev. A 30, 1185 (1984). 

[9] L.D. Landau and E.M. Lifshitz, Electrodynamics of con- 
tinuous media (Pergamon Press, Oxford, 1960). 
[10] L. Onsager, Phys. Rev. 37, 405 (1931); 38, 2265 (1931); 



9 



[11] W. Eckhardt, Opt. Commun. 41, 305 (1982); Phys. Rev. 

A 29, 1991 (1984). 
[12] W.C. Chew, Waves and Fields in Inhomogeneous Media 

(IEEE Press, New York, 1990). 
[13] J.D. Jackson, Classical Electrodynamics, 2ndedn. (Wiley, 

New York, 1975). 
[14] E.M. Purcell, Phys. Rev. 69, 681 (1946). 
[15] C.-T. Tai, Dyadic Green Functions in Electromagnetic 

Theory (IEEE Press, New York, 1993). 



[16] Z. Xiang and Y. Lu, IEEE Trans. Microwave Theory 
Tech. 44, 614 (1996). 

As pointed out in Ref . 0| , this paper contains some crit- 
ical mistakes. The relation between the scattering coeffi- 
cients and the Green tensor in this paper is inconsistent. 

[17] L.-W. Li, M.-S. Leong, T.-S. Yeo, and P.-S. Kooi, J. of 
Electromagnetic Waves and Appl. 14, 961 (2001). 

[18] In Ref. Q, the angular factor is erroneously given as 1/10. 



